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Abstract. A short overview is given on the development of our present paradigm of the large scale structure of the Universe 
with emphasis on the role of Ya. B. Zeldovich. Next we use the Sloan Digital Sky Survey data and show that the distribution of 
phases of density waves of various scale in the present-day Universe are correlated. Using numerical simulations of structure 
evolution we show that the skeleton of the cosmic web was present already in an early stage of the evolution of structure. The 
positions of maxima and minima of density waves (their phases) are the more stable, the larger is the wavelength. The birth 
of the first generation of stars occured most probably in the central regions of rich proto-superclusters where the density was 
highest in the early Universe. 
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INTRODUCTION: YA. B. ZELDOVICH 
AND THE DEVELOPMENT OF OUR 

UNDERSTANDING OF THE 
LARGE-SCALE STRUCTURE OF THE 
UNIVERSE 

This conference is devoted to Ya. B. Zeldovich, to cele- 
brate his 95th birthday, so let me start with an overview 
of our contacts with Ya. B. Zeldovich and his group, in 
particular of his role in the development of the modern 
view of the large scale structure and the evolution of the 
Universe. Thereafter I will analyse the luminosity den- 
sity field of galaxies of the Sloan Digital Sky Survey 
(SDSS). Finally I will use numerical simulations of evo- 
lution of structure to follow the changes of phases (posi- 
tions of maxima and minima) of density waves of various 
scales, using wavelet techniques. 

Tartu astronomers met Ya. B. Zeldovich for the first 
time in 1962. The main building of the new observa- 
tory in Toravere was not yet completed, but under the 
initiative of Ya. B. Zeldovich a summer school on cos- 
mology was held there, one of the first of similar sum- 
mer and winter schools. Fig. 1 shows participants of 
the school visiting the old Tartu Observatory. In 1960s 
Tartu astronomers participated in Astrophysics Seminars 
in Sternberg Institute, organized by Ya. B. Zeldovich, 
during their visits to Moscow. Here astronomers from all 
Moscow astronomical institutes gathered to discuss ac- 
tual problems of astrophysics. 

My own first report in this seminar was in 1971, when 
I made a summary of my doctoral thesis on galactic mod- 
eling. This report was a starting point of our collabora- 
tion with Zeldovich and his group. Next year members 
of our cosmology group were invited to participate in 
the Archoz Winter School on astrophysics, where I made 



a more detailed report of new galactic models. In Jan- 
uary 1974 the School was held near the Elbrus moun- 
tain in a winter resort. I reported our very fresh results 
on the study of the dynamics of satellite galaxies. Our 
data suggested that all giant galaxies have massive in- 
visible coronas, which exceed masses of known stellar 
populations at least tenfold. Thus dark matter must be 
the dominating component in the whole universe. In the 
discussion after the talk two questions dominated: What 
is the physical nature of the dark matter? and What is its 
role in the evolution of the Universe? 

The importance of dark matter for cosmological stud- 
ies was evident, thus Tartu astronomers organized in Jan- 
uary 1975 a conference in Tallinn devoted solely to dark 
matter (Doroshkevich et al. [9]). Historically this was the 
first conference on dark matter. In addition to dynamical 
and morphological properties of galaxy systems (Einasto 
et al. [13, 14]) the nature of galactic coronas was dis- 
cussed. Here are titles of some talks held at this con- 
ference: Zeldovich [45]: "Deuterium nucleosynthesis in 
the hot Universe and the density of matter'' \ Komberg 
& Novikov [24]: ''Physical nature of galactic coronas'' \ 
and Jaaniste & Saar [21]: ''Properties of stellar halos" . 
These studies reported at the conference demonstrated 
that both the stellar as well as gaseous models for coro- 
nas have difl&culties. It is very difl&cult to explain physical 
properties of the stellar corona, since all old spheroidal 
populations have a low mass-to-luminosity ratio whereas 
this ratio for the corona is very high. Also, no fast- 
moving stars as possible candidates for stellar coronas 
were found. Gaseous corona cannot consist of neutral gas 
(intergalactic hot gas would ionize the coronal gas), but 
a corona consisting of ionized gas would be observable. 
Actually, at least part of the corona must be gaseous, 
as demonstrated by the morphology of satellite galax- 
ies: elliptical (non-gaseous) companions lie close to the 



FIGURE 1. Participants of the Cosmology Summer School 1962 visiting old observatory in Tartu. 



primary galaxy whereas spiral and irregular (gaseous) 
companions of the same luminosity have larger distances 
from the primary galaxy, as shown by Einasto et al. [14]. 
The amount of the hot gas is, however, insufficient to ex- 
plain flat rotation curves of galaxies. 

Whatever the nature of galactic coronas, it was evident 
that dark matter as the dominant population of the uni- 
verse determines details of physical processes of the for- 
mation and evolution of galaxies. Zeldovich asked Tartu 
astronomers for help in solving the question: Can we find 
observational evidence that can be used to discriminate 
between various theories of galaxy formation? At this 
time there were two basic rivaling theories of structure 
formation: the pancake theory by Zeldovich [44], and 
the hierarchical clustering theory by Peebles [30] (see 
Figures 2 and 3). According to the Zeldovich scenario 
the structure forms top-down: first matter collects into 
pancakes and then fragments to form smaller units. In 
the hierarchical clustering scenario the order of the for- 
mation of systems is opposite: first small-scale systems 
(star-cluster sized objects) form, and by clustering sys- 
tems of larger size (galaxies, clusters of galaxies) form; 



this was called the bottom-up scenario. 

When solving the Zeldovich question we recalled our 
previous experience in the study of galactic populations: 
structural properties of populations hold the memory of 
their previous evolution and formation. Random veloci- 
ties of galaxies are of the order of several hundred km/s, 
thus during the whole lifetime of the Universe galaxies 
have moved from their place of origin only by about 
1 Mpc (we use in this paper the Hubble constant 
in the units of Hq = 100 h km s~^ Mpc~^). In other 
words - if there exist some regularities in the distribu- 
tion of galaxies, then these regularities must reflect the 
conditions in the Universe during the formation of galax- 
ies. What is expected in the pancake model is shown in 
Fig. 4. This simulation of structure formation was made 
in mid-1970s in the Institute of Applied Mathematics 
with the largest available computer at this time by the 
Zeldovich group. In the simulation a network of high- 
and low-density regions was seen: high-density regions 
form cells which surround large under-dense regions. 

When we started this study, first all- sky complete red- 
shift surveys of galaxies were just available: the Shapley- 




FIGURE 2. Zeldovich with his wife visiting Estonia. 



Adams revised catalogue (Sandage & Tammann [32]), 
complete up to the magnitude 13.5, and the Second 
Revised Catalogue of Galaxies (de Vaucouleurs et al. 
[8]), containing data on even fainter galaxies. The three- 
dimensional distribution of galaxies, groups and clusters 
of galaxies can be visualised using wedge diagrams. In 
these diagrams, where galaxies as well as groups and 
clusters of galaxies were plotted - see Fig. 5, regularity 
was clearly seen: galaxies and clusters are concentrated 
to identical essentially one -dimensional systems, and the 
space between these systems is practically empty. 

The distribution was quite similar to the distribution 
of simulation particles in a numerical simulation of the 
evolution of the structure of the Universe prepared by 
the Zeldovich group. Fig. 4. We had the impression that 
observed high-density regions - superclusters of galaxies 
- could be identified with Zeldovich pancakes. 

We reported our results (Joeveer & Einasto [20]) at 
the lAU symposium on Large-Scale Structure of the Uni- 
verse in Tallinn 1977, the first conference on this topic. 
The main results were: (1) galaxies, groups and clus- 
ters of galaxies are not randomly distributed but form 
chains, converging in superclusters; (2) the space be- 
tween galaxy chains contains almost no galaxies and 
forms holes (voids) of diameters of 20 . . .50 Mpc, as 
seen from the wedge diagrams in Fig. 5. The presence 
of holes (voids) in the distribution of galaxies was re- 
ported also by other groups: Tully & Fisher [38], Tiff't 
& Gregory [36] and Tarenghi et al. [35] in the Local, 
Coma and Hercules superclusters, respectively. Theoret- 
ical interpretation of the cellular structure was discussed 
by Zeldovich [46] . 

The most important aspect of the structure of super- 
clusters is the fine structure of chains: clusters and groups 




FIGURE 3. Jim Peebles and Scott Tremaine at the lAU 
Symposium on Large Scale Structure of the Universe in Tallinn 
1977. 



FIGURE 4. Distribution of particles in a N-body simulation 
by the Zeldovich group in 1975. 



within a chain are elongated along the chain, as seen in 
the long chain of clusters, groups and galaxies of the 
Perseus-Pisces supercluster. This chain is located almost 
perpendicular to the line of sight, and the scatters of po- 
sitions of clusters/groups along the chain in radial (red- 
shift) and tangential directions are practically identical. 
This demonstrates that the chain is essentially an one- 
dimensional structure. A direct consequence of this ob- 
servation is that galaxies and groups/clusters of the chain 
are already formed within the chain. A later inflow from 
random locations to the chain is excluded, since in this 
case it would be impossible to stop galaxies and clus- 
ters in the chain after the inflow. This is possible only 
in a gaseous dissipative medium before the formation of 





FIGURE 5. Wedge diagrams for two declination zones. Filled circles show rich clusters of galaxies, open circles - groups, dots 
- galaxies, crosses - Markarian galaxies. In the 15° - 30° zone two rich clusters at RA about 12 h are the main clusters of the 
Coma supercluster, in the 30° - 45° zone the clusters at RA about 3 h belong to the main chain of clusters and galaxies of the 
Perseus-Pisces supercluster. Note the complete absence of galaxies in front of the Perseus-Pisces supercluster, and galaxy chains 
leading from the Local supercluster towards the Coma supercluster in the upper zone and towards the Perseus-Pisces supercluster 
in the lower zone (Joeveer & Einasto [20]). 



galaxies. 

New data gave strong support to the Zeldovich pan- 
cake scenario. However, some important differences be- 
tween the model and observations were evident. First of 
all, there exists a rarefied population of simulation par- 
ticles in voids, absent in real data. This was the first in- 
dication for the presence of physical biasing in galaxy 
formation - there must be primordial gas and dark mat- 
ter in voids (because gravity cannot evacuate voids com- 
pletely), but due to low matter density no galaxy forma- 
tion takes place here. 

The second diff'erence lies in the structure of galaxy 
systems in high-density regions: in the original pancake 
model large-scale structures (superclusters) have rather 
diff'use forms, real superclusters consist of multiple in- 
tertwined filaments (Joeveer & Einasto [20], Zeldovich 
et al. [47], Oort [27]). The reason of this discrepancy was 
clarified several years later. 

Let me recall one more historical moment in our 
search for understanding the cosmic evolution in late 
1970s. At the Tallinn symposium Parijskij [29] made a 
report on his search for primordial perturbations in the 
Universe with RATAN-600, the largest and most sensi- 
tive radio telescope at this time. No temperature fluc- 
tuations were found, the upper limit was about 10""^ of 
the mean temperature. After the talk Zeldovich discussed 
these results with Parijskij and expressed his doubts that 



something must be wrong in his observations. Theoret- 
ical calculations show that at the epoch of recombina- 
tion the density (and temperature) fluctuations must have 
an amplitude of the order of 10"^, otherwise the struc- 
ture cannot form, since gravitational instability that is re- 
sponsible for the growth of the amplitude of fluctuations, 
works very slowly in an expanding Universe. 

The above arguments concern fluctuations of the bary- 
onic gas density. Then astronomers recalled the possible 
existence of non-baryonic particles, such as heavy neu- 
trinos. This suggestion was made independently by sev- 
eral astronomers (Cowsik & McClelland [7], Szalay & 
Marx [34], Tremaine & Gunn [37], Doroshkevich et al. 
[11], Chernin [6], Bond et al. [5]) and others. They found 
that if dark matter consists of heavy neutrinos, then this 
helps to explain the paradox of small temperature fluc- 
tuations of the cosmic microwave background radiation. 
This problem was discussed at a conference in Tallinn in 
April 1981. Recent experiments by a Moscow physicist 
Lyubimov were announced, which suggested that neu- 
trinos have masses. If so, then the growth of perturba- 
tions in a neutrino-dominated medium can start much 
earlier than in a baryonic medium, and at the time of 
recombination perturbations may have amplitudes large 
enough for structure formation. At the conference ban- 
quet Zeldovich gave an enthusiastic speech: ''Observers 
work hard in sleepless nights to collect data; theorists 



interpret observations, are often in error, correct their 
errors and try again; and there are only very rare mo- 
ments of clarification. Today it is one of such rare mo- 
ments when we have a holy feeling of understanding the 
secrets of Nature." Non-baryonic dark matter is needed 
to start structure formation early enough. 

However, in the neutrino-dominated dark matter 
small-scale perturbations are damped. This prevents the 
formation of fine structure, such as galaxy-sized halos 
and filaments. Motivated by the observational problems 
with neutrino dark matter, a new dark matter scenario 
was suggested by Blumenthal et al. [1], Bond et al. 
[4], Pagels & Primack [28], Peebles [31], Bond & Sza- 
lay [3], Doroshkevich & Khlopov [10] with hypothetical 
particles as axions, gravitinos, photinos or unstable neu- 
trinos playing the role of dark matter. This model was 
called the Cold Dark Matter (CDM) model, in contrast 
to the neutrino-based Hot Dark Matter model. 

Numerical simulations of the evolution of structure for 
the hot and cold dark matter were compared by Melott 
et al. [26], and by White et al. [40, 41] (the standard 
model with the density parameter = 1, SCDM). In 
contrast to the HDM model, in the SCDM scenario for- 
mation of structure starts early and superclusters consist 
of a network of small galaxy filaments, similar to the ob- 
served distribution of galaxies. Thus SCDM simulations 
reproduce quite well the observed structure with clusters, 
filaments and voids, including its quantitative character- 
istics (percolation or connectivity, the multiplicity distri- 
bution of systems of galaxies). 

As a further step to improve the model, the cosmologi- 
cal constant. A, was incorporated into the scheme. Argu- 
ments favouring a model with the cosmological constant 
were suggested already by Gunn & Tinsley [19], Turner 
et al. [39], Kofman & Starobinskii [23]: combined con- 
straints on the density of the Universe, ages of galaxies, 
and baryon nucleosynthesis. Numerical models with the 
cosmological A-term were developed by Gramann [18], 
a post-graduate student of Enn Saar. Comparison of the 
SCDM and ACDM models shows that the structure of 
the cosmic web is similar in both models. However, in 
order to get a correct amplitude of density fluctuations, 
the evolution of the SCDM model has to be stopped at 
an earlier epoch, thus the ACDM model is superior of 
the SCDM model. The ACDM model combines essen- 
tial aspects of both original structure formation models, 
the pancake and the hierarchical clustering scenario. 

One difficulty of the original pancake scenario by 
Zeldovich is the shape of objects formed during the 
collapse. It was assumed that forming systems are flat 
pancake-like objects, whereas dominant features of the 
cosmic web are filaments. This discrepancy has found 
an explanation by Bond et al. [2]. They showed that in 
most cases just essentially one-dimensional structures, 
i.e. filaments form. Two coauthors of this study, Kofman 



and Pogosyan, are graduates of the Tartu University and 
former members of our cosmology group. 



DISTRIBUTION OF GALAXIES IN THE 
SLOAN SURVEY 

Presently the largest project to map the Universe is the 
Sloan Digital Sky Survey (SDSS) (York et al. [42], 
Stoughton et al. [33], Zehavi et al. [43]). The goal is to 
map a quarter of the entire sky: to determine the positions 
and photometric data in 5 spectral bands for galaxies and 
quasars (about 100 million objects) down to the red mag- 
nitude r = 23, and the redshifts of all galaxies down to 
r = 17.7 (about 1 million galaxies). All seven data re- 
leases have been made public. This has allowed to map 
the largest volume of the Universe so far. 
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FIGURE 6. The luminosity density field in spherical shells 
of 10 Mpc thickness and distances of 240 Mpc and 
325 Mpc is shown in the upper and lower panels, respec- 
tively. In the upper panel the very rich supercluster in the lower 
region is SCL126 (Einasto et al. [15]), or the Sloan Great Wall. 
Cellular arrangement of superclusters and filaments surround- 
ing voids is clearly seen. The angular diameters of cells in the 
325 /^-i Mpc shell are smaller than in 240 h ^ Mpc shell. 

The luminosity density fields of two spherical shells 
of the full Sloan Survey in the Northern hemisphere are 
shown in Fig. 6. The luminosity densities were corrected 
for the incompleteness due to the distance dependent 
lower luminosity limit of the flux-limited survey. The 
cosmic web with filamentary superclusters and voids 
is very well seen. Faint filaments crossing large voids 



can be seen also. The web has a characteristic cellular 
pattern, the diameter of one cell - a void surrounded 
by superclusters - is fairly constant, about 100 Mpc. 
Since shells are located at different distances from us, the 
angular diameters of faraway cells are smaller. The great 
variability of supercluster richness can be easily seen. 

The variability of supercluster richness is due to joint 
action of density perturbations of different scales. To un- 
derstand better the role of perturbations of various scale 
we have decomposed the density field into components 
of various scale using a wavelet technique. 

As a first step we analysed a thin, almost 2- 
dimensional wedge of the Sloan survey equatorial slice, 
2.5 degrees thick, centered on the rich supercluster 
SCL126 in the catalog by Einasto et al. [15] (the 
Sloan Great Wall). The high-resolution luminosity den- 
sity field, reduced to a constant thickness by dividing the 
local density by the relative thickness of the wedge at a 
particular distance, is shown in the upper left panel of 
Fig. 7. 

The importance of phase information in the formation 
of the supercluster-void network has been understood 
long ago, as demonstrated by Szalay (personal commu- 
nication) by randomizing phases of a Voronoi simulated 
picture of a cellular network. To study the role of phase 
information in more detail we Fourier transformed the 
density field and randomized phases of all Fourier com- 
ponents and thereafter Fourier transformed it back to see 
the resulting density field which has for all waves the 
same amplitudes (power spectrum) as the original field; 
see the upper right panel of Fig. 7. We see that the whole 
structure of superclusters, filaments and voids has gone, 
the field is fully covered by tiny randomly spaced density 
enhancements - groups of galaxies, there are even no rich 
clusters of galaxies in this picture that were comparable 
in luminosity to the real clusters. 

Next we used wavelet analysis to investigate the role 
of density waves of different scale (for details of the 
wavelet analysis see Martinez & Saar [25]). The full field 
was divided into 8 components, so that the component 
with the largest waves represents waves of the length 
about 512 Mpc (the size of the box), and each of the 
following components represents waves of twice smaller 
scale. Actually, each component contains waves of scale 
interval from ^f2 times lower to V2 times higher than 
the mean wavelength, so that by combining waves of 
all components we restore accurately the initial density 
field. Fig. 7 shows components 4 to 7; they characterize 
waves of length from about 32 to 256 Mpc. 

Let us discuss the distribution of waves in more detail. 
The middle left panel of Fig. 7 shows the waves of length 
about 256 Mpc. The regions of highest density of this 
component are three very rich superclusters: located in 
the upper left part of the Figure, near the center, and 
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FIGURE 7. The upper left panel shows a rectangular region 
extracted from the high-resolution luminosity density field of 
the Northern equatorial slice of SDSS, and converted to dis- 
tances using the dimensionless Hubble constant h = 0.8. The 
density was calculated using Gaussian smoothing with a rms 
scale of 0.8 Mpc. The observer is located at the lower left 
corner. The upper right panel shows the same density field as 
in the previous panel, but with randomly shifted phases. The 
following panels show the wavelet analysis of the slice. The 
middle left panel shows density waves of scales approximately 
equal to half of the box size, each following panel shows waves 
of scale 2 times smaller. 



in the lower right part of the panel, the superclusters 
SCL155, SCL126 and SCL82 in the Hst by Einasto et al. 
[15], respectively. The middle right panel shows waves 
of scale of about 128 Mpc. Here the most prominent 
features are the superclusters SCL126 and SCL82 as in 
the previous panel, but we see several new superclusters: 
N15 and N04 in the Hst by Einasto et al. [12]. The 
supercluster Nl 5 lies between SCL126 and SCL155 near 
the minimum of the density wave of 256 Mpc scale. 

The lower left panel plots waves of scale about 64 
Mpc. Here all superclusters that are seen at larger scales 
are also visible. A large fraction of density enhancements 
are either situated just in the middle of the low-density 
regions of the previous panel, or they divide massive su- 
perclusters into smaller subunits. The highest peaks are 
substructures of massive superclusters, and there are nu- 
merous smaller density enhancements between the peaks 
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FIGURE 8. The high-resolution density field of the model Ml 00. Left panel shows the field at the redshift z = 100, the right 
panel at the present epoch (z = 0). The color coding is chosen to show regions of overdensity, i.e. regions of local density greater 
than unity in mean density units. 



of the previous panel. 

This analysis allows to make several conclusions. First 
of all, our results show that in all cases superclusters form 
only in regions where large density waves combine in 
similar phases of high-density regions to generate high 
density peaks. Very rich superclusters are objects where 
density waves of all large scales (up to the wavelength 
- 250 Mpc) have similar phases of density maxima. 
The smaller the maximum wavelength of such phase syn- 
chronization, the lower is the richness of superclusters. 
Similarly, large voids are caused by large-scale density 
perturbations of wavelength ~ 100 Mpc, here large- 
wavelength modes combine in similar phases of minima 
of density waves to generate under- densities. 

The second conclusion that can be made from the 
wavelet analysis concerns the relation of waves of dif- 
ferent scale. Density perturbations which have random 
phases and wavelengths in respect to the principal pertur- 
bations (responsible to the formation of rich and very rich 
superclusters) cancel partly each other out. An excep- 
tion are the perturbations which have similar phases and 
wavenumbers that are integer multiples of the wavenum- 
ber of the principal perturbation. These are the harmon- 
ics (overtones) of the principal perturbation. The first 
overtone has a wavenumber twice of the principal one; 
it forms a density peak just in the middle of the void 
between two rich superclusters, forming two sub- voids 
of approximately equal diameter within the large under- 



dense region. This behavior is well seen in the wavelet 
decomposition of the density field in Fig. 7. 

The phase synchronization that generates rich super- 
clusters is not surprising, if considered from the point of 
view of the density field. A high peak in the field is pos- 
sible only in the case if density waves of various scale 
up to a certain limiting wavelength have maxima at the 
same location. The acoustic behavior of density waves - 
the domination of overtones - deserves more attention. 
Presently it is not clear how to explain this behavior. 

THE EVOLUTION OF 
SUPERCLUSTERS ON SIMULATIONS 

To investigate the role of density perturbations 
of diflTerent wavelength we used N-body simu- 
lations with 256^ particles and with box sizes 
L = 64, 100, 256, 768 h'^Mpc. Most models 
were calculated by Ivan Suhhonenko in Tartu, the 
L=100h-^ Mpc model was obtained by Volker Miiller 
in Potsdam; we denote this model as Ml 00. In the 
present review I shall analyse only the model Ml 00, 
other models yield similar results. In model Ml 00 
the following cosmological parameters were used: the 
matter density = 0.27, the baryonic matter density 
= 0.04, the dark energy density ^Ia = 0.73; the initial 
amplitude parameter erg = 0.84, and the dimension- 




less Hubble parameter h = 0.72. Calculations started 
at an early epoch, z = 100, the early evolution was 
calculated using the Zeldovich approximation. Particle 
positions and velocities were extracted for the redshifts 
z = 100, 10, 5, 3, 2, 1, 0.0. For every particle we 
calculated also the local density in units of the mean 
density, using positions of 27 nearby particles. 

For all models and epochs we calculated the high- 
resolution density field using the B3 spline smoothing 
with the step width equal to the spatial resolution of the 
simulation (for details of smoothing with the B3 spline 
see Einasto et al. [16]). The high-resolution density field 
of the Ml 00 model at the epochs z = 100 and z = 
is presented in Fig. 8. Here we have selected from the 
full simulation box a sheet that passes through a large 
void surrounded by superclusters and filaments. At the 
high-redshift epoch the density contrast is much smaller 
as expected, but main features of the filament system 
around the void are clearly seen. The basic high-density 
knots (protohalos of present-day rich clusters) are very 
well seen already at the early epoch, and their positions 
have changed only very little. 

The Figure 9 shows the wavelet decomposition of 
the same field. Left panels show the high-redshift, and 
right panels the present epoch decomposition. Wavelets 
of the order from 6 (upper panels) to 4 (lower panels) 
are shown, they show waves from the size of the box 
(100 Mpc) to about 25 Mpc. Only overdensity 
regions are shown, having densities above unity in mean 
density units. The colour code is chosen so that it corre- 
sponds approximately to the factor of linear growth from 
z= 100 to z = 0, about 76. 

The Figure shows that large-scale waves have almost 
the same shape in the early and present epoch, only the 
amplitude has grown as expected. Even the waves of 
scale 25 Mpc have changed only little - the main 
features at both epochs are the same, only some filaments 
have joined during the evolution. This is the well-known 
eff'ect of structure merging inside superclusters. 

It is well known that the Zeldovich approximation of 
linear growth is the more accurate the earlier epochs and 
larger waves we consider. Thus it is evident that at epochs 
earlier than our first epoch discussed here, the growth can 
be described very accurately by the Zeldovich approx- 
imation. In other words: all large-scale structures have 
had almost no evolution in their positions of density wave 
features (maxima, minima and other similar features), 
the evolution of these features is only in the amplitude 
of the density contrast. 

The second important lesson from this analysis con- 
cerns the positions of maxima of waves of various scale. 
When we take a density maximum of the wave of the 
largest scale, then at the same position ALL waves of 
smaller scale have also a maximum, i.e. maxima of 
density waves are synchronized. The same behavior is 




FIGURE 9. Wavelet analysis of the model Ml 00 shown in 
Fig. 8. Left panels show density waves of various scales at the 
redshift z = 100, right panels at the present epoch (z = 0). In 
the upper panels the characteristic wavelength is approximately 
equal to the size of the box, L = 100 Mpc, in the middle 
panels the scales are twice smalle, and in the lower panels four 
times smaller. 



present in the real density field of Sloan galaxies as dis- 
cussed above. 

Thus wavelet analysis of our simulations demonstrates 
that positions of high-density objects (superclusters, rich 
clusters) almost do not change during the whole history 
of the evolution of structure. In other words: the skeleton 
of the cosmic web or the supercluster-void network was 
fixed already at a very early stage, much before the epoch 
of recombination. This conclusion confirms early ana- 
lytic calculations by Kofman & Shandarin [22], where 
main contours of the cosmic web were found using the 
theory of adhesion at the epoch which corresponds to the 
end of the inflation stage and to the start of the Hubble 
expansion. 



Our analysis shows that high-density peaks of the den- 
sity field are present already at the very early stage of the 
evolution of the structure. Recent numerical studies show 
that the formation of first generation stars starts just in- 
side these high-density peaks. The density is highest in 
protohalos of central clusters of future rich superclusters. 
Simulations by Gao et al. [17] have shown that metal- 
free gas in dark matter haloes of virial temperature about 
2000 K and mass M ~ 10^ M© cools efficiently, thus giv- 
ing rise to formation of stars. In these high-density re- 
gions star formation can start as early as at the redshift 
z ~ 50. These first generation stars have large masses and 
evolve rapidly, and at least some of them may explode 
and spread products of nuclear synthesis to the surround- 
ing "pure" gas. 

CONCLUSIONS 

• Superclusters form in regions where large-scale 
density waves combine in similar locations of max- 
ima. 

• Superclusters are the richer the larger is the wave- 
length of wave synchronization. 

• Medium rich superclusters are often formed near 
minima of waves of scale ~ 250 Mpc by the first 
overtone of this wave. 

• Voids form in regions where large-scale density 
waves combine in similar locations of minima. 

• Large-scale density waves do not change positions 
of maxima during the evolution, i.e the web of the 
filaments that defines the final state of the large- 
scale galaxy distribution is actually present in the 
initial density inhomogeneities. 

• In centers of rich superclusters all waves of various 
scale have maxima. Since positions of large waves 
do not change, high local densities were formed 
already at a high redshift z> 100. Thus formation 
of the first generation stars started in the centers of 
rich superclusters. 
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